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Abstract 

We introduce the notion of Differential Sequences of ordinary differential equations. 
This is motivated by related studies based on evolution partial differential equations. 
We discuss the Riccati Sequence in terms of symmetry analysis, singularity analysis 
and identification of the complete symmetry group for each member of the Sequence. 
The singularity analysis reveals considerable structure for the values of the coefficients 
of the leading-order terms and resonances of the different principal branches. Full 
proofs of the symmetry properties are performed for differential equations defined by 
their recursion properties and not given in explicit form. 



1 Introduction 



In the study and application of differential equations there are certain equations 
which are prominent by means of their utility of application or their ubiquity of occur- 
rence. In the particular case of ordinary differential equations examples of members 
of this prominent class are the Riccati equation, the linear second-order equation, the 
linear third-order equation of maximal symmetry, the Ermakov-Pinney equation, the 
Kummer-Schwarz equation and its generalisation. These equations are not indepen- 
dent. As Conte observed [BJ, the study of any of the Riccati equation, the linear 
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second-order equation and the Kummer- Schwartz equation is equivalent to a study of 
the other two. More recently the closely connected linear third-order of maximal sym- 
metry and the Ermakov-Pinney equation have been added to the list given by Conte. 
In [lUj the generalised Kummer-Schwartz equation has also been included. Since these 
equations are of different orders and have differing symmetry properties, it is evident 
that the connections among them are nonlocal. 

A considerable part of the motivation for this study is found in the works of Pe- 
terssen et al [20] and Euler et al [8] in which the authors obtain recursion operators for 
linearisable 1 + 1 evolution equations. The specific equation of relevance to this work 
is the eighth of their classification, videlicet 



where As is an arbitrary constant and h 8 is an arbitrary function of the dependent 
variable u. The corresponding recursion operator is 



The usual area of application of recursion operators has been that of partial differ- 
ential equations and in particular evolution equations of which the equation above is a 
sample. Euler et al [9] applied the idea of a recursion operator to ordinary differential 
equations. They did not present a treatment of all possible classes, but concentrated 
upon two representative equations, videlicet the Riccati equation and the Ermakov- 
Pinney equation. 

The recursion operator which generates the Riccati Sequenced is 



where D denotes total differentiation with respect to x. The recursion operator es- 
sentially originates from Case VIII, R$[u], for partial differential equations when the 

2 In [5] the term 'Riccati Hierarchy' is used. For a reason which becomes apparent in what follows we 
believe that the term 'sequence' is more appropriate. The 'correct' term to be used to describe this sequence 
has been a matter of perception of the language. 



Ut = u xx + X 8 u x + h$u 2 x 



VIII 



R 8 [u] = D x + h 8 u x . 



R.O. = D + y 
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dependence on t is removed, ie (ut = 0), and one writes u x = y. We note that in 
applications to ordinary differential equations the notation usually used in the context 
of partial differential equations can be and is simplified. 

When (pQ) acts orjfl y (which can be denoted as Ro) it generates, in succession, all 
the members of the Sequence which we denote by R n , videlicet 

y' + y 2 = (2) 

y" + 3yy' + y 3 = (i? 2 ) (3) 

y'" + Ayy" + 3y' 2 + 6y V + y 4 = (i? 3 ) (4) 

y"" + 5yy"' + Wy'y" + Wy 2 y" + I5yy' 2 + Wy 3 y' + y 5 = (i? 4 ) (5) 

(D + y) n I/ = (i? n ). (6) 

As an aid in our discussion we distinguish between an element of the Sequence, denoted 
as indicated above by R n , and its left side by writing the latter as R n . 



In addition to the formal definition given in (|6|) there is a differential recurrence 
relation. 



Lemma: The members of the Riccati Sequence satisfy the differential recurrence rela- 
tion 

d 

— (Rrn+i) = (m + 2)R m . (7) 



Proof: By inspection of ([2]), (J3]) and (j4|) it is evident that ([7]) is true for the initial 

3 This is a matter of notation. In this draft we have chosen the simplest. In another context, especially 
when one wishes to be 'physically' correct, one would choose cxp[J ydx] since, if y is replaced by —x one 
obtains the generating function for the solution of the time-independent Schrodinger equation for the simple 
harmonic oscillator. 
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members of the sequence. We assume that 

d 

— (R k+1 ) = (k + 2)R k (8) 

for some k > 2. Then 

(D + y)R k+ i 

= [(D + y)d y + 1] Rk+i 

= (D + y)(k + 2)R k + R k+l 

= {k + 3)R k+1 

as required. 

Remark: This recurrence relation is a little intriguing bearing in mind the connection 
to the generating function for the Hermite polynomials. One recalls the differential 
recurrence relation for the Hermite polynomials, videlicet 

dH n (x) 

— = 2nH n -x(x), 

ax 

and notes the similarity of structure to that given in ([7|). 




The Riccati Sequence contains as its two lower members the Riccati equation and 
the Painleve-Ince equation. The Riccati equation [22] has an history now of almost 
three centuries and is one of the few examples known for which there exists a nonlinear 
superposition principle. Its close relationship to the linear second-order differential 
equation and the Kummer-Schwartz equation has already been noted [6]. The Painleve- 
Ince equation [191 [T2] is a notable example of a nonlinear second-order differential 
equation of maximal symmetry possessing both Left and Right Painleve Series 
and arises in a remarkable number of applications^- 



In this paper we report the notable properties of the members of the Riccati Se- 
quence. We concentrate upon the number of Lie point symmetries, singularity prop- 
erties, first integrals, explicit integrability and complete symmetry groups. In Section 

4 For a listing see [171 [TT1 [5l [3], The authors of the third paper place (O in a broader class of equations 
which they describe as being of modified Emden type. 
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2 we present the symmetry analysis of the members of the Sequence. In terms of 
the explicit integrability of the members of the Sequence the lack of point symmetry 
in the higher-order members is remarkable and indicates the necessity for considering 
nonlocal symmetries. We note that this was the case in the treatment of a pair of 
equations of Ermakov-Pinney type [15]. In Section 3 we elaborate upon the results 
of the singularity analysis of the lower members of the Riccati Sequence presented by 
Euler et al [9] by considering the pattern of the values of the resonances for the general 
member of the Sequence. Section 4 is devoted to a consideration of the invariants and 
first integrals of the members of the Sequence. A short discussion on the symmetries 
of the integrals is presented. In Section 5 we present the complete symmetry group of 
the general member of the Sequence after firstly considering the results for to give 
a concrete basis for the theoretical discussion. In Section 6 we recall the solution of 
R n and extend our discussion to an equation containing combinations of the R{. In 
the case of the second-order equation so formed the properties have been known for a 
long time \19\ \12\ 113] . Among our concluding remarks in Section 6 we summarise the 
remarkable properties found for the Riccati Sequence which is based upon a rather el- 
ementary recursion operator. We indicate that the route to complexification presented 
in Section 6 broadens the class of differential equations which may be subsumed into 
a more general concept of a Riccati Sequence. 



2 Symmetry analysis 



As a first-order ordinary differential equation ([2]) possesses an infinite number of 
Lie point symmetries. Equation ([3]) was examined in [T7] for its Lie point symmetries 
which were found to be the following eight, indeed the maximal number for a second- 
order ordinary differential equation, here written in a more elegant way to match the 
results of the succeeding members of the sequence, 

ri = x 2 yd x - y[2 + xy{xy - 2)]3 y 
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T 2 = yd x - y 3 d y 

T 3 = xyd x + y 2 (l - xy)d y 

T 4 = xd x - yd y 

T 5 = x 3 {xy-2)d x -x{xy-2)[2 + xy(xy-2)}d y 

T 6 = -x 2 (xy - 2)d x + xy(xy - 2){xy - l)d y 

T 7 = x 2 d x + (2 - 2xy)d y 

r 8 = d x 
with the algebra sl(3,R). 

We proceed with R3 to find an unexpected three, ie 

Ti = d x 

T 2 = xd x - ydy 

T 3 = x 2 8 x + (3 - 2xy)d y 

with the algebra sl(2,R). 

Proposition: The general member of the Riccati Sequence possesses the symme- 
tries 

Ti = d x 

T 2 = xd x - ydy 

T 3 = x 2 d x + (n- 2xy)d y 

with the algebra sl(2,R). 

Remark: For R2 there are an additional five symmetries and for Ri an additional 
infinity. 

Proof: The possession of Ti is straightforward since ([6|) is autonomous. We proceed 
with r 2 . The nth extension of T 2 is 

n 

T^=xd x -^2(j + l)y^d yU) (9) 

j=0 
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and it is a simple calculation to demonstrate that 



-2R X 
-3R 2 - 



We assume that the statement 



(10) 



is correct and we prove that 



rlf^-Rm+l — —{m + 2)Rm+l- 



We commence with the statement 

T^R m+1 = T^ +1 \D + y)R ri 



(11) 



We break the calculation of the right side of (jlip into some parts. Firstly we have 



r [m+l] ;jD ' 



LB 



+ DT [ l n+l] . 



(12) 



The Lie Bracket is 

m+l 

xd x + l)y®d y v) A + E y {k+1) d yW 

j=0 fc=0 
'"m+l 



LB 



E(i + 1 )^' } ^)'Ey (fc+1) V 

j=0 k=0 
m+l 



- dx -Y^(j + i)Y,[y {3)d y^^ k+1)d : 

j=0 k=0 
m+l 

A-E^' + 1 )E{y ^ 



(k) 



LB 



LB 



(fc) - y 



j=0 fc=0 
m+l 



-d x ~ U + 1) {v U) dyU-V ~ y {j+1) 9yU) } 



3=0 
m 



m+l 



-dx - E0' + % a+1) ^> + E V + i)y°' +1) ^) 

j'=o i=o 
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3=0 

= -D + (m + 2)y( m+2 ~>d y{m+lh 

where the overdot on the summation means that the sum is taken to whatever order 
of derivative required. 

In a much simpler calculation we have 

T^ +1] yR m = {T^ +1] y )R m + y (T^ +1] R m ) 

= - y R m + y rl m+1] R m . (13) 

We substitute the results of both the Lie Bracket and (113j) into (lllj) . We have 
if +11 i? m+ i = \-D + (m + 2)y^d y{m+1) + DT [ ™ +1] 



= -(D + y)R m + (D + y)T^ +1] R 

® -(D + y^-p + ^(771 + 1)^ 

= -(m + 2)(D + y) J R m 

= -(m + 2)^+1 (14) 



7/1 



which proves the result. 



Note that in the proof we have made use of the fact that R m contains derivatives 
only up to y( m \ 

We turn now to the proof in the case of 1^ . 

It is a simple matter to show that 

if 1 ^ = -2x(l + l)Ri 
Yfk 2 = -2x(2 + l)i? 2 
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rf ] i? 3 = -2x(3 + l)R 3 
and one can easily assume that 

T [ ™ ] R m = -2x(m + l)R m 
so that the task is now to demonstrate that from this property it follows that 

4 m+11 R m+1 = -2x(m + 2)R m+1 . 



(15) 



(16) 



We commence with the definition and so we have 

ri m+1] iW = T [ ™ +l \D + y)R m 

Tf +1 \{D + y)\ R m + (D + y)T [ ™ +1] R m 

J LB 

rf (D + y)} R m + {D + y) (rf 1 + d y ) R m , 

where the last line is a consequence of remembering the definition of I3 and that the 
(m + l)th derivative in if +1 1 does not act on R m . 

The first task is to compute the Lie Bracket. We have 

m+l m+l 



fc=0 



% + (m + 1 - 2xy)d y -2xY, U + l)2/ 0) ^o) - ]T j(j + l)y {j ~ 1] d yU) , 

3=1 3=1 
m+l 

-2xd x + 2yd y + 2 Y, (J + l)y (j) d y m + 2xy'd y + (m + 1 - 2xy) 

3=1 

m+l 

+2x U + 1) E [v {k+1)d y^ > V (j)d yt 
j=l k=0 
m+l 



3=1 



k=0 
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We compute the two subsidiary Lie Brackets separately. 

m+l 

LBi = 2x (j + 1) E {v lk+1)s kjd y v) - y®8 jM id v m } 



m+l 



fc=0 



2x £(j + 1) {y^fi^, - 



2x < 



(m + 2)y( m+2 )^ ( „ 1+1) - 2y'd y - £ y^"^ 



W) 



In a similar way we find that 

m 

LB 2 = (m + l)(m + 2)y^d y{m+1) - 2yd y -2^(j + l)y®d yU) . 

3=1 

We return to the main calculation and after some simplification of terms which cancel 
we continue as below. 



p[m+l]^ 

3 



m+l 



-2x\d x + y'd y + Y,y {j+1) d l 

5=1 



U) 



2(m + 2)/ m+1 )^ (m+1) 



+m + 1 - 2xy + (m + l)(m + 2) 2 / m+1 ^ Cm+ i ) 
+2x(m + 2)y( m+ Vd y(m+1) + (£> + y) {-2x(m + 1) + fy)} i? m 
{-2x(D + y) + (m + l) + (D + y)d y - 2(m + 1)(-D + y)x} 
-2x(m + 2)(D + y)£ m + \{D + - (m + 1)] i? m 
-2x(m + 2).R m+1 



since 



[(£> + yH - (m + 1)] R m = 0. 

This last statement is equivalent to the result of Lemma [H ie 

d 

-Q- (Rm+i) = (m + 2)R m 

and so the result is proven. 

Remark: Note that in the proof above we have again made use of the fact that R n 
contains derivatives only up to y( m \ 
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3 Singularity analysis 



We use singularity analysis as a tool to determine whether a given differential 
equation is integrable in terms of functions almost everywhere analytic. For equations 
which pass the Painleve Test we are then encouraged to seek closed- form solutions. 
We find that all members of the Riccati Sequence possess the Painleve Property. As 
we see below, explicit integrability follows. Independently of this integrability the 
results of the singularity analysis provide some very interesting patterns in terms of 
the parameters, that is the p, a and r of the standard analysis. Similarly interesting 
patterns have been reported in [16j . 

Euler et al [9] presented the singularity analysis of the Riccati Sequence. For all 
elements the leading-order behaviour is ax -1 with the possible values of a being listed 
in the Tables. We summarise their results. 

Table I be here 

In Table II we advance from [9] and present the properties in terms of the singularity 
analysis for the general member of the Riccati Sequence. In that way we are able to 
comment upon the integrability or otherwise of all members of the Sequence in the 
sense of Painleve. 

Table II be here 

The pattern of the resonances is as follows: The set of resonances for a = j is 
obtained from the set for a = j — 1 by subtracting the number n + 1 from the largest 
positive resonance of the latter set. 
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It follows from Tables I and II that all members of the Riccati Sequence pass the 
Painleve Test and therefore each is integrable in terms of analytic functions. 



When we apply the Riccati transformation 



V = a- (17) 

w 



on the nth member of the Riccati Sequence, we observe that the most simplified equa- 
tion, ie ■u/ n+1 ) = 0, arises when we chose a = 1 which is a consequence of the singularity 
analysis itself [3]. Therefore (fTT|) may be written as 



x = x, w = exp 



/ 



ydx 



It is a matter of simple calculation to verify the following proposition. 

Proposition: The general solution of the nth member of the Riccati Sequence, 
n > 1, is given by 

Vn- v n a i , 

where the Ai, i = 0,n, are constants of integration. 



4 Invariants and First Integrals 



For the purposes of this paper we commence with some definitions. 

Definition 1: An invariant I of an ordinary differential equation is said to be any 
nontrivial function which satisfies 

on solution curves of the differential equation. 
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Definition 2: A first integral of an nth-order ordinary differential equation is any 
function / = I(x, y,y' , . . . , which satisfies (|19p . 

Note that we do depart from some standard definitions for reasons which become 
evident below. There are various conventions concerning the meanings of the expres- 
sions 'first integral', 'invariant' and 'conserved quantity'. We are not concerned by the 
third in this paper since we do not use that expression. Some writers use the first two 
expressions interchangeably. Others prefer to distinguish between the two by insist- 
ing that the former be autonomous whereas the latter is allowed to depend upon the 
independent variable. This is a sensible distinction under appropriate circumstances. 
Indeed in terms of the integration of an ordinary differential equation the distinction 
can be quite critical. In this paper we have varied the definitions to suit the very pre- 
cise purpose of distinguishing between two classes of function both of which have the 
property of having a zero total derivative on solution curves of the differential equation. 

We wish to identify all invariants and first integrals, as defined above, of each 
member of the Riccati Sequence. In order to do so we take the rath member of that 
Sequence, perform an increase of order by using the Riccati transformation to obtain 
w (n+l) = o -^y e com p U te the fundamental integrals and invariants of that equation 
and by reverting the transformation one can deduce the whole set of invariants for the 
rath member of the Riccati Sequence, videlicet 

( j \ 
7 J = ( E J-— J xJ ~ iR (n-i) J exp 

where = 1 by convention. 
Note that for the invariants of w^ n ' = we would have the above invariants (|20p except 
that instead of R( n -i) we would have w^ n+1 ~ 1 ^ and there would be no exponential. 

We now turn our attention to first integrals. One takes the ratio of two separate 
invariants to obtain a first integral. An independent set for R n can be variously defined. 



ydx 



1,71+ 1, 



(20) 
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The set {Tij}, defined by 

Fij = y, j = 1,» + l,n, (21) 

-n 

is an independent set of first integrals for R n . Such a simple formula is not available if 
one wishes to describe a set of autonomous first integrals. 

An interesting aspect arises in the symmetry properties of the first integrals which 
we briefly note. If one computes the symmetries of all first integrals of ([3]), one finds 
that they all share the algrebra A\ © s A2- By a curious misfortune this feature does 
not persist and therefore one is confronted with an absolute zero for contact (not to 
mention point) symmetries for the first integrals of 



5 Complete symmetry groups 



The concept of a complete symmetry group of a differential equation was introduced 
by Krause [H] as the group associated with the set of symmetries, be they point, 
contact, generalised or nonlocal, required to specify the equation or system completely. 

We start with to give a flavor of the procedure and then we prove the general 
result for any member of the Riccati Sequence. 



Proposition: The complete symmetry group of is given by the symmetries 



A.1 



A, 



exp 



exp 



ydx 



ydx 



{y}dy 
{xy - l}dy 



5 When equation (j4|) is related to = by an increase of order, the symmetries of the first integrals 
of the latter become, by inverting the transformation, nonlocal symmetries for the ones of the former. 
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A 3 
A 4 



exp 
exp 
exp 



ydx 
ydx 
ydx 



{x 2 y - 2x}d y 
{x 3 y - 3x 2 }d y 
{x A y - 4x 3 }dy. 



Proof: We calculate the fourth extensions of Ai — A5. 



A 



exp 



ydx 



{yd y + (y' - y 2 )d y , + {y" - 3yy' + y 3 )d y 



+(y"' - ±yy" - 3y' 2 + 6yV - y A )d y ,„ 

+(y {lv) - 5yy m - Wy'y" + Wy 2 y" + 15yy' 2 - 10y 3 y' + y 5 )d y(lv) } 



A 



exp 



ydx 



{-d y + 2yd y , + 3{y' - y 2 )d y „ + A(y" - 3yy' + y 3 )d y 



+5(y"' - Ayy" - 3y' 2 + 6y 2 y' - y 4 )cW + xA^ 



A 



-2 exp 



ydx 



{-dy, + 3yd y „ + Q(y' - y 2 )d y i, + 10(y" - 3yy' + y 3 )d y(lv) ) 



+2xA^ ff + x 2 A^ 



AM 



-6 exp 



J ydx {-d y » + Ayd y w + W(y' - y 2 )d y{iv) } 



+3^ + 3*^ + *^ 



—24 exp — J ydx 
+4xA' 4 J // + 6x 2 a' 4 J„ 



{-dyin + 5yd {i „)} 



3e// 



+ . 4 a[ 4 ], 



where the subscript ejff stands for the effective part of that symmetry. 

When we act all the above extensions on the general fourth-order ordinary differential 

equation, videlicet 



y 



we obtain the system of five equations 



W=f(x,y,y>,y"y>), 



iiv) 



Of 



df 



5yy"' - Wy'y" + Wy 2 y" + Ibyy' 2 - Wy 3 y' + y^ = y^- + (y' - y 2 )^ 



dy' 
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+(/ - 3yy' + y 3 )|4 + (</" ~ %/ - 3y' 2 + 6yV - z/)#4 (22) 



dy'' 



5(y"> - 4yy" - 3y' 2 + 6yV - y 4 ) = + 2y^ 
W-y 2 )|4+4(/-3yy / + y 3 )-' 9 



10(y" - 3W + y 3 ) = -f£ + 3y^ + 6(y' -y 2 ) 9 



dy" 
dy' ' " a dy 1 ' 



dy'' 



5y 



which can be solved backwardsly to give the following expressions for the derivatives 
of all arguments in / 



df_ 
dy'" 

df_ 

dy" 
dj_ 

dy' 
df 
dy 



When these expressions are substituted into <\22h we recover R^. 



-by 

-10(y' + y 2 ) 

-10(2/" + Syy' + y 3 ) 

-5(y'" + 4yy" + 3y /2 + 6yy + y 4 ). 



Proposition: The complete symmetry group of R n is given by the (n + 1) symme- 



tric 



A; 



exp 



ydx 



{x^y -{i- l)x l ~ 2 }d y , i = 1, n + 1 



Proof: We remind the reader that R n stands for the nth member of the Riccati 
Sequence and R n for the left hand side of R n . It is essential for the notation required in 



6 The correct number of symmetries required to specify an equation completely has been discussed 
and 0. 
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this proof to introduce the adjoint recursion operator, R.O. a = D — y, which generates 
the adjoint Riccati Sequence, videlicet 

y'-y 2 = (R?) (23) 

y"-3yy' + y 3 = (R%) (24) 

y'»-4yy"-3y' 2 + 6y 2 y'-y 4 = (R%) (25) 

yW - 5yy"' - Wy'y" + lOy 2 y" + Ihyy' 2 - Wy V + y 5 = (R a 4 ) (26) 
y W _ QyyM _ I5y' y '" + i^y'" + 60 y? /y" - lOy" 2 - 2(W' 

+15y /3 - 45y 2 y' 2 + 15yV - 2/ 6 = (27) 

p-yfy = (K)- (28) 
The left hand side of the above members is denoted by R®- 



Remark: We note that there is a reflection here between the Riccati Sequence and 
its adjoint in that the formulas for the Sequence are reflected in the formulas for the 
adjoint Sequence. 

Lemma: The general members of the Riccati and the adjoint Riccati Sequences satisfy 
the recurrence relation 

Rn = R° + ^ ( . jRf^Rn-i- (29) 
i=l ^ ' 



Proof: It is easy to show that (|29p holds for n = 2 and n = 3. We assume that ([29 
is true for n = k and we prove that it is true for n = k + 1. 



R 



k+l 



(D + y)R k = (D + y) 
k 



i=l 



Rt +l + 2yR% + J2[ J { D ft-i**-* 



R-t + ^l + ; jR?-\Rk-i 

i=l ^ ' 

yRf-iRk-i^ 



k + l 
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i=i ^ 1 ' 

+ib( k+1 ) (^-i+i - ^-o + 2/ e ( k + x ) ^-i^-* 

i=l ^ % ' i=l ^ 1 ' 



{n = ] R a k+l + 2yR% + iz{ k + l ) k ? k ^ + it( k + 1 ) RtiRk-i + i + y(Rk - R\ 

i=l ^ 1 ' i=l ^ 1 ' 

= + y-Rfe + y-Rfe + E L _ i ) ^f-i^-i+i + E 

i=2 ^ Z ' i=l 



fc + 1 



-Rf_i-Rfc-i+i 



fe+i 



Rt +l+ yR% + yR k + Y J \ (^l) + (T 1 ) 



i=l 



fe+1 



i=l 



which proves the result. 



The nth extension of A,- is 



A 



exp 



ydx 



(i - i)!0 w (i- a) + ^ydyii-D + {l ^-R?a y{i 



2! 



(^+2)!~ (n + 1)! ~ Q 

H K 2 v (i+i) + ■■■ + — „■ , N, -Kn-i+lP y (n) 



3! 

2 1 d fc 



fc=0 



where 



A (l)e// = A l '^(-i) = 0. #„(<>) = fy) ^(i) = etc > 



— — — etc 

dx° ' dx 1 da; 



The action of the effective part of A , 



n—i+l 



A L 1 



-a y( i-2) + ^ 



j=0 



on the general nth-order ordinary differential equation, videlicet 

y^=f(x,y,y',y",y"',...,y^), 
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gives 



n + 1 



ji-i + 2 
For i = 1 (l30l) gives 



n— i+l 



9/ + 8/ 



j=0 



n-1 



i = 1, n 



j=0 



Of 
dyV) ' 



The remaining n of (|30p are equivalent to 



n + 1 U Q df 
n-i + l) n ~ l 



+ E 



5/ 



■j— i— 1 



* = 1, 



System (j32[) may be conveniently represented as the system of equations 



where 



(F„)i 
(Qn)y 



n+1 

i 

g/ 
.7 



i = l,n, 



We define 



(Qn) 



and 



(L r 



n+1 



, i = l,n, 

Rj-i-l, i<3, 

-1, * = j, 

0, i>j. 

Rj-i-i, % < j, 

-1, i = j, 

0, i>j. 

Rn—ii % — 1, 71. 



We prove the following two Lemmas. 
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Lemma: The matrices and Q n defined in fa34\ ) and A35\) are the inverses of each 
other, ie they satisfy 

QnQn = I„- (36) 



Proof: That ([36]) is true for n = 1 is obvious. We assume that it is true for n = k, ie 



QZQk = i k . 



(37) 



Then 



Qfc+iQfc+i 







Qfe — Lfc 


or -i. 




o£ -i 


QfcQfc, - 








1 



ED 



provided Q?L/% + L? = 0. We consider the ith element of this term. 



J k)i 



(QgL fc + L^ = ^(Q^-(L jfc ) J - + (L 



fc + 1 

i 



Rk~i-Rk-i- XT ( + _ • ) R( j-i-l R k-3 



0. 



Lemma: 



QnL^ 



(38) 
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Proof: Consider an element of the product Q„L°. This is 

0' 



(Qn.)i?(L")j 



n + 1 

i 



Rn-i + 

j=i+l 

n+1' 



Rj-i 



■j-i-l 



J 



R 



n-j 



n+1 

i 



n+1 



Rn-i + y] 

j=i+l 



x R n _jRj—i—i 



j\ fn + 1\ /n + 1 



.7 



(L n )i. 



Therefore (1331) is written as 



F n — (Qn) 

BP n T» 

— MA 

^ -L n . (39) 

Equation (j3Tj) is equally written as 

9/ 



pa \ r £?« 



which, when (J59"|) and (|29j) are used, results to (i2 n ) and the proposition is proven. 



6 Concluding remarks 

We are familiar with the use of sequences of numbers and functions in the math- 
ematical and wider scientific literature. The last several decades have seen intensive 
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study of hierarchies based upon nonlinear evolution equations which have their bases 
in the mathematical modelling of physical phenomena. The extension to ordinary dif- 
ferential equations has only been made recently [9]. We have chosen to term these 
related equations as differential sequences to reflect the dual nature of a definition 
which combines the idea of an operation which generates the elements of the sequence 
and also that the elements of the sequence are composed of derivatives. It seems to us 
that the word sequence is more appropriate in this context than hierarchy and is more 
in keeping with the tradition of mathematical terminology. 

The Riccati Sequence, which has been the subject of the present study, illustrates 
with a degree of excellence, which one hopes can be surpassed, the type of mathematical 
properties which are likely to make such sequences an object of fond study. The 
Riccati equation has a long and distinguished history in both the theory of differential 
equations and the application to divers phenomena. The Riccati Sequence and its 
adjoint are based on recursion operators closely related in form to the Dirac operators 
of the quantum mechanical simple harmonic oscillator which in themselves are simply 
autonomous versions of two of the Lie point symmetries of the classical simple harmonic 
oscillator. With such a lavish heritage one should not be too surprised that the Riccati 
Sequencj^l exhibits such properties in terms of symmetry and integrability. 

Each element of the sequence can be linearised by means of the nonlocal transfor- 
mation - often called a Riccati transformation - and so is trivially integrable. Although 
Ri and R2 - the Riccati and Painleve-Ince equations - display exceptional symmetry in 
the sense of Lie point symmetries, the remaining elements of the sequence possess just 
the three-element algebra sl(2, R) of Lie point symmetries. The distinguishing feature 
of the Lie symmetries of the elements of the sequence is the possession of n + 1 (in the 
case of R n ) exponential nonlocal symmetries which completely specify the elements. 
The algebra of these symmetries is (n + 1)A±. 



7 Equally the adjoint sequence. In the following discussion the properties of the adjoint sequence may be 
inferred mutatis mutandis from those of the Riccati Sequence. 
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As a consequence of the ability to linearise the equations of the Riccati Sequence 
each element possesses an invariant derived from the so-called fundamental integrals 
of the parent linear equation^. These invariants are not integrals in the conventional 
sense as they contain the integral of the dependent variable. However, the integral 
of the dependent variable appears as an exponential term and so functions of these 
invariants which are homogeneous of degree zero in the dependent variable are first 
integrals. It was for this reason that we introduced specialised meanings for the two 
expressions, 'invariant' and 'first integral', for the purposes of this paper. Since there 
are (n + 1) linearly independent (exponential nonlocal) invariants, there are n func- 
tionally independent integrals which reflect most adequately the integrability of each 
member of the Riccati Sequence. 



In terms of singularity analysis the Riccati Sequence possesses properties which 
may even be regarded to outshine the symmetry properties. Not only is each element 
of the Sequence explicitly integrable in terms of analytic functions apart from isolated 
polelike singularities - a property which has been used to illustrate certain subtle and 
not widely appreciated implications of singularity analysis [3] - but also there are 
patterns to the possible values of the coefficients of the leading-order terms and the 
values of the resonances for each of the principal branches. In passing one recalls that 
each branch is principal and the need to be concerned with the implications of branches 
with unfortunate properties is obviated. The patterns of the resonances deserve further, 
separate treatment, particularly in respect of other possible sequences. 



The Riccati Sequence has provided an excellent vehicle for the introduction of 
the notion of differential sequences of ordinary differential equations. Its Recursion 
Operator is simple. Its generating function is elementary. Two of the elements, Ri 
and i?2; are well-known in the literature as well as in applications. Already Euler et al 
[9] have made a brief excursion into the properties of the Ermakov-Pinney Sequence, 



8 Thcse integrals have been termed fundamental since an ?2th-ordcr linear ordinary differential equation 
possesses n linearly independent integrals which are linear in the dependent variable and its first (n — 1) 
derivatives. The integrals are known for their interesting algebraic properties. 
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as the name suggests, based upon the Ermakov-Pinney equation [21] which has been 
found in many varied applications. 

In the Introduction we have recalled the classical association of the Riccati equation, 
the linear second-order ordinary differential equation and the third-order Kummer- 
Schwarz equation. Since these equations are related by means of nonlocal transforma- 
tions, there is no reason to exclude from this select group other equations which are 
simMy rel a te cg A.ea* in W are initial studies of a sequence based upon the 
Kummer-Schwarz equation [TU] for which the Recursion Operator is now an integrodif- 
ferential operator of third order by way of contrast with the simple Recursion Operator 
of the Riccati Sequence. One looks forward to further revelations of fascinating prop- 
erties of other differential sequences. 

As a final remark we recall that the evolution equations and associated Recursion 
Operators [20] which are the source material for this work contain arbitrary parameters 
and unspecified functions. In this work we have deliberately kept to the minimum 
required to make a sensible sequence. It is a little like treating the autonomous linear 
oscillator rather than a time-dependent, damped and forced oscillator. The latter can 
be transformed to the former by means of well-defined transformations and a study 
of the former is easier due to the simplicity possible in the presentation. Despite our 
personal preference for simplicity we do nevertheless accept that there are those who 
wish to see the treatment of the more general systems. 



Consider 

n 

E n = J2fi( t ) R i = °- (40) 

8=0 

Equation (|40p is linearisable by means of the Riccati transformation (|17|) . with a = 1, 
to 

n 

X>(i)^ +1) =0. (41) 

i=0 



9 One is well aware that all properties may not travel well through the process of nonlocal transformation, 
but that obstacle is already encountered under far less exotic transformations. 
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Equation (|4ip possesses (n + 1) solution symmetries and the homogeneity symmetry. 

Whilst there can be no dispute that the solution of (|4ip is more difficult than that 
of u/ n+1 ) = 0, for our purposes there is no difference. Subject to some mild conditions 
on the functions fi(t) - basically that they be continuous [12] [p 45] - the solutions of 
(|4ip exist and that is all that is required for inclusion into the general framework of 
the treatment above. An nth-order linear ordinary differential equation can have n + 1, 
n + 2 or n + 4 Lie point symmetries which represent the n solution symmetries, linearity, 
autonomy (in the right variables) and the general sl(2,R) symmetry of equations of 
maximal symmetry [18]. There is no gainsaying that the additional symmetries make 
the process of solution in closed form just that much easier. However, all we need is 
the existence of the n solution symmetries to provide the symmetries necessary for our 
purposes. Consequently we have taken the clearer part so that the ideas and results 
contained in this paper be as evident as can be possible in such matters. 

We look forward to a development of the subject of differential sequences of ordinary 
differential equations to match that of the related areas in partial differential equations. 
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Table 1: Singularity analysis for the first four members of the Riccati Sequence. 



Member 


Leading-order coefficients 


Resonances 


Ri 


a = 1 


r = — 1 


R2 


a = 1 
a = 2 


r = -1,1 
r = -1,-2 


Rs 


a = 1 
a = 2 
a = 3 


r = -l,l,2 
r = -l,l,-2 
r = -1,-2,-3 


i?4 


a = 1 
a = 2 
a = 3 
a = 4 


r = -l, 1,2,3 
r = -1,1,2,-2 
r = -l,l,-2,-3 
r = -1,-2, -3, -4 
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Table 2: Singularity analysis for the general member of the Riccati Sequence. 



Member 


Leading-order coefficients 


Resonances 




a = 1 


r = -l,l,2,...,n-l 




a = 2 


r = -l,l,...,n-2,-2 




a = n 


r = —1, —2, . . . , — n 
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